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Fermion bags, duality and the three dimensional massless lattice Thirring model 
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The recently proposed fermion bag approach is a powerful technique to solve some four-fermion 
lattice field theories. Due to the existence of a duality between strong and weak couplings, the 
approach leads to efficient Monte Carlo algorithms in both these limits. The new method allows us 
for the first time to accurately compute quantities close to the quantum critical point in the three 
dimensional lattice Thirring model with massless fermions on large lattices. The critical exponents 
at the quantum critical point are found tohe v = 0.85(1), = 0.65(1) and r;^ = 0.37(1). 
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Strongly interacting quantum critical points contain- 
ing massless fermionic excitations are interesting from 
many perspectives. Such critical points describe second 
order phase transitions in strongly correlated two dimen- 
sional materials such as Graphene [ij. The search for 
critical points in gauge theories in three spatial dimen- 
sions containing massless fermions is interesting in the 
context of the physics beyond the standard model [1]. 
Properties of these fermionic quantum critical points re- 
main poorly understood due to lack of reliable compu- 
tational techniques. While the renormalization group 
approach based on large N or e-expansion techniques 
can provide a glimpse of the rich possibilities, accurate 
Monte Carlo (MC) calculations are desirable. Unfortu- 
nately, MC methods for lattice field theories with mass- 
less fermions in three or more space-time dimensions con- 
tinue to be challenging. For example, the popular Hybrid 
Monte Carlo (HMC) method [3] either suffers from sign 
problems or encounters severe singularities due to small 
eigenvalues of the fermion matrix. All HMC calculations 
performed so far have always relied on extrapolations to 
the massless limit which are known to be difficult [3| . As 
far as we know MC calculations close to a quantum criti- 
cal point with exactly massless fermions on large lattices 
do not exist. 

Recently, a new approach called the fermion bag ap- 
proach was proposed as an alternative method to solve 
a class of lattice field theories with exactly massless 
fermions . It is an extension of the meron cluster idea 
proposed some time ago @ . The idea behind the fermion 
bag is to identify fermion degrees of freedom that cause 
sign problems and collect them in a bag and sum only 
over them. This is in contrast to traditional approaches 
where all fermion degrees of freedom in the entire thermo- 
dynamic volume are summed to solve the sign problem. 
When the fermion bag contains only a small fraction of all 
the degrees of freedom and the summation can be per- 
formed quickly, the fermion bag approach can be used 
to design powerful MC methods. Sometimes, the bag 
splits into many disconnected pieces further simplifying 



the calculation. 

The general idea of the fermion bag can be illus- 
trated easily with the following example. Consider lat- 
tice fermion models formulated with 2n Grassmann vari- 
ables per lattice site denoted as iJi{x) and ^^{x), where 
i = 1, 2, ..n represent flavor indices and x denotes the Eu- 
clidean space-time lattice point containing V sites. Let D 
be the V xV free fermion matrix whose matrix elements 
are denoted as D^y- We will assume that the proper- 
ties of D are such that the following fc-point correlation 
function involving the flavor i : 

Ci{xi, ...,Xk) = J [dV#]exp(^^ i^iix) D^y ipi{y) 

tpi{xi)ipi{xi) ... ■^p^{xk)■^p^{xk) (1) 

is always positive. An example of such a matrix D, is 
the massless staggered fermion Dirac operator which is 
popular in constructing four-fermion lattice models 
It is easy to prove that 



a{xi,..,Xk)=Bet{D) Det(G[{x},]) 



(2) 



where G[{a;}i] is the k x k matrix of propagators be- 
tween the k sites Xp,p — 1 , . . , fc whose matrix elements 
are Gx^^x, = D~^^ . It is also possible to argue that [sj. 



C,{xi, ..,Xk) = ^ct{W[{x},]) 



(3) 



where the matrix W^[{a;}i] is a (V^ — fc) x (V" — fc) matrix 
which is the same as the matrix D except that the sites 
{x} = {xp,p = l,2,..fc} are dropped from the matrix. 
The identity 



T)ei{D) Det(G[{a;}i]) = Det(W[{a;},]) 



(4) 



leads to a concept of duality in the fermion bag approach 
as we explain below. 

Consider a generic four-fermion lattice field theory ac- 



strong coupling fermion bags 



Weak coupling fermion bag 



FIG. 1. An illustration of a "fermion-bag" configuration at strong couplings (left) and weak couplings (right). The interactions 
in this illustration are represented by solid bonds and the fermion bags are represented by shaded region. At strong couplings 
the fermion bag is made up of free sites and breaks up into many disconnected pieces, while at weak couplings the bag contains 
interaction sites. 



tion involving n massless staggered fermions given by 

5 = - ^ D^yij,{y) 

x,y,i 
{xy} ij 

where (xy) refers to some well defined set of neighboring 
lattice sites. Further we will assume that all the couplings 
Uij^(xy) ^I'c non- negative real constants. Many interest- 
ing four- fermion models are of this type [7J . The partition 
function can be expanded in powers of the coupling and 
is given by 

/n 
e-^ = ^{[C/f} l[C,ixi,...,Xk,), (6) 
[{^}] «=i 

where {[?7]}^ refers to a generic power of the coupling 
and ki^i = 1,2, ..n refers to the number of interaction 
vertices for each flavor. On a finite lattice the expansion 
is convergent since it is a polynomial. 

The above expansion of the partition function begs 
for the following intuitive interpretation. Since the ki 
interaction sites contain both tpi and -011 the j"^ fla- 
vor of fermions are already paired on these sites and do 
not cause sign problems. On the other hand, unpaired 
fermions of the i*'^ flavor that hop freely on the remain- 
ing sites can indeed cause sign problems and need to be 
summed over to solve the sign problem. The free sites 
are collectively referred to as a fermion bag. The sum- 
mation of all fermion world lines inside the bag leads to 
the weight Ci{xi,X2, --^Xk) = Det{W[{x}i]) which is the 
determinant of a (V^ — ki) x (y — ki) matrix. This deter- 
minant can be evaluated easily if (V — ki) is small. This is 
expected at strong couplings and hence we refer to these 



free fermion bags as strong coupling fermion bags. It was 
shown in Q that at strong couplings a fermion bag splits 
into many small disconnected pieces making things even 
simpler. The left figure of Fig. [1] gives an illustration of 
the disconnected pieces of a strong coupling fermion bag. 

At weak couplings the above definition of a fermion 
bag loses its charm since V — ki becomes large. How- 
ever, thanks to a concept of duality we can construct 
the fermion bag differently. At weak couplings we can 
view the interactions as the unpaired fermionic degrees 
of freedom that cause fluctuations over the paired free 
fermionic vacuum. In this view the fermions hop from 
one interaction site to another interaction site leading 
to sign problems and need to be summed over. Now 
the interaction sites form the fermion bag. Again the 
summation of the fermions inside this dual bag leads 
to the same weight Ci{xi,X2, --fXic) = Det(M^[{x}i]) = 
Det(i?) Det(G[{a:}i]) but now viewed as the determinant 
of a fci X ki matrix. Note we have used the duality relation, 
Eq.Q here. The determinant can now be calculated eas- 
ily since ki is small at weak couplings. Hence we refer to 
these dual bags as weak coupling fermion bags. The right 
figure of Fig. [1] gives an illustration of a weak coupling 
fermion bag. The weak coupling fermion bag approach is 
equivalent to the idea of summing over all Feynman dia- 
grams and was introduced earlier in the framework of di- 
agrammatic determinantal Monte Carlo method Q . On 
the other hand, in our opinion the fermion bag approach 
is more intuitively appealing in the context of lattice field 
theories since it uncovers the powerful concept of duality 
and extends to strong couplings. 

The fermion bag approach is general and can 
be adapted to relativistic Wilson fermions and non- 
relativistic fermions. However, in some models the weight 
of a fermion bag is no longer a determinant, but involves 
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new mathematical structures like fermionants [9|, which 
can be exponentially difficult to compute [lo| . In such 
cases the fermion bag approach is not useful. However, 
these models can be changed by introducing unconven- 
tional interactions (like six or higher fermion interac- 
tions), while still preserving many interestingsymmetries 
and the corresponding low energy physics [5|- In these 
exotic models, the fermion bag weight is again a deter- 
minant and the fermion bag approach becomes useful. 
Thus, many new fermion models can be solved with the 
fermion bag approach. 

As a first application of the new approach we have 
studied the three dimensional massless lattice Thirring 
model with two Grassmann valued fields per site denoted 
by ip{x) and il^{x). The action is given by 

ip{x) D^y ip{y) tp{x)il;{x) ■4;{y)ip{y), 



x,y 



{xy) 



(7) 

where D is the massless staggered fermion matrix, {xy) 
refers to the set of nearest neighbor sites of a cubic lattice 
and U is the coupling that generates the current-current 
coupling of the continuum Thirring model. We use 
anti-periodic boundary conditions in all three directions. 
The lattice model is invariant under a Uf{l) x U^{\) 
symmetry, where Uf{l) is the fermion number symme- 
try and U^{1) is the chiral symmetry. When the cou- 
pling is small the model contains four flavors of mass- 
less two-component Dirac fermions at long distances due 
to fermion doubling. When the coupling is large, chiral 
symmetry breaks spontaneously and generates a single 
massless Goldstone boson and the fermions become mas- 
sive. Thus, the model contains a quantum critical point 
Uc which separates a phase with massless fermions from a 
phase with massless bosons. The quantum critical point 
has been studied earlier using mean field techniques [lH , 
and traditional MC methods [l2l-[l5jl. A variant of our 
lattice model has been used recently to study a related 
quantum critical point in the context of Graphene 16| . 

Close to the quantum critical point a continuum quan- 
tum field theory description of the long distance physics 
must emerge. This continuum theory should contain four 
flavors of two component Dirac fermions in three Eu- 
clidean dimensions. As was discussed in [l^, the lat- 
tice interactions generate many continuum four-fermion 
interaction terms and the continuum Lagrange density 
takes the form 

^ = V'i(a;)(CT • V)V'i(a;) + \^gA\i^i(x)a{TA)ijil^j{x)\'^ 
+ ~gA\^,{x){TA\j^l^,{x)Y] (8) 

where ipi{x),^pi{x),i = 1, ..,4 are the four flavors of two 
component Dirac fermion flelds, a are the three Pauli 
matrices, Ta,A = 1,..,16 are the sixteen generators of 
the U{4:) group in the flavor space under which the free 
theory is invariant. Repeated indices are assumed to be 



summed over. The couplings gA and gA are only con- 
strained by the lattice symmetries. They take values 
such that the 17(4) symmetry of the free theory is broken 
to a Uf{l) X U^{1) subgroup [l3|- As far as we know, 
a renormalization group (RG) flow analysis in this rela- 
tively large but constrained space of couplings, within a 
controlled approximation such as large N or e-expansion, 
is not available and should be an interesting topic for fu- 
ture research. However, the existence of a quantum crit- 
ical point in the lattice model does imply that the RG 
analysis will find a nontrivial fixed point with one rele- 
vant direction. Here we compute the critical exponents 
at this fixed point through MC calculations. 

The fermion bag approach for the above lattice model 
was first developed in [5| and it was shown that the par- 
tition function can be written as 



Z = ^U^''r>et{W[n]) 



(9) 



where [n] refers to the configuration of Nh interaction 
bonds and W[n] is the {V - 2Nb) x {V - 2Nb) staggered 
Dirac matrix restricted to the free sites. These free sites 
form the strong coupling fermion bag (see left figure in 
Fig. [1]). At the quantum critical point Ni, is about an 
eighth of the lattice volume and hence the above form 
of the partition function is not useful. However, thanks 
to duality we can think in terms of the weak coupling 
fermion bags (see right figure of Fig. [Ij. For the above 
model, the duality relation (Eq. (|3])) takes the form 



(10) 



where G[n\ is an iVf, x Ni, free fermion propagator ma- 
trix between even and odd lattice sites of the interaction 
bonds. Using recent algorithmic advances fl8|. we have 
constructed an efficient determinantal Monte Carlo algo- 
rithm for this problem [Toj . If the autocorrelation times 
and equilibration times are measured in units of a sweep, 
we flnd that our algorithm has no further critical slowing 
down even at the quantum critical point. 

In order to uncover the properties of the quantum crit- 
ical point wc focus on three observables (Let L be the 
lattice size): The chiral condensate susceptibility (20j . 



the chiral winding number susceptibility 



(11) 



(12) 



(defined through the conserved chiral charge {qx)a — 
J2x£S ^x Vx,a {D~^)x,x+a + Z^xss ^Ea; passiug through 
the surface S perpendicular to the direction a, where 
£x = (1)^1+^2-1-3:3 _ ^]^g staggered fermion phase factors 



4 




FIG. 2. Plots of (^x) a-iid RfL'^^'^'^ as a function of [/ for L from 12 to 40. The solid lines show the combined fit 

which give Uc — 0.2608(2), i/ — 0.85(1), rj — 0.65(1) and rj^ = 0.37(1). The inset shows a linear-log plot of Rf versus L at U = 
(solid line) and U — 0.260 (dashed line). The lines are fits to a power law. 



Vx,a = exp(z7r^a 
Ci = (0,0,0), 6 
ratio 



x) are defined through the 3-vectors 
= (1,0,0) and 6 = (1,1,0)), and the 



i?/ =Cf(L/2-l)/Cf(l), 



(13) 



(defined through the fermion two point function Cp (d) = 
^J2'L=i{'^3:'^x+d&) in which x belongs to a site with 
Ex — and d is a unit vector along each of the three 
directions). Since the fermions are exactly massless, in 
the vicinity of Uc we expect these three observables to 
satisfy the following simple finite size scaling relations: 



k=0 
3 



k=0 
3 



R 



(14a) 
(14b) 
(14c) 



fc=0 



where we have kept the first four terms in the Taylor 
series of the corresponding analytic functions. Our goal 
is to compute the critical exponents 77, v and 77^ at the 
quantum critical point. Large N calculations usually give 
V = rj — \ and 77^ = Q. 

Earlier studies of the quantum critical point were fo- 
cused on computing v and 77. In these studies the four 
fermion coupling is converted into a fermion bilinear us- 
ing an auxiliary field. The fermions are then integrated 
out and the remaining problem is solved using the HMC 
method. In order to avoid singularities, all calculations 
are done with a finite fermion mass. A close examination 
of the earlier work reveals that different analysis produce 
substantially different results. However, the presence of 
large errors makes everything look consistent. In our 
opinion, the presence of two infrared scales (the fermion 
mass and the length of the box) makes the analysis diffi- 
cult. In contrast, since we work with massless fermions, a 



single combined fit to Eqs. p4)) with sixteen parameters 
is easy. Indeed, a combined fit of all our data from 12"^ to 
40^ lattice gives us = 0.85(1), 77 = 0.65(1), 77^ = 0.37(1) 
and Uc = 0.2608(2) with a x^/d.o.f. = 1.3. The complete 
list of the sixteen fit parameters are listed in Table HI Our 
computation of 77^ is new. 



/o 



h 



/2 



h 



2.52(3) -2.53(5) 0.71(3) 0.10(1) 0.369(3) 0.63(1) 



K2 



PO 



Pi 



P2 



PZ 



0.52(2) 0.09(1) 33.9(2) -5.0(1) -2.0(2) -2.5(5) 

TABLE I. Results from the combined fit of the data to 
Eqs. p4|l . In addition to the above twelve parameters, the 
fit also gives Uc = 0.2608(2), ?? = 0.65(1), v = 0.85(1), 
= 0.37(1) with a xVd-o.f of 1.3. 



Plots of our data along with the fits are shown in Fig. [5] 
Since our data fits very well to the expected scaling form 
for a whole range of lattice sizes, we feel confident that 
the corrections to scaling are small. 
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